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Abstract 

This study addresses the issue of oil removal from water using hydrophilic porous membranes. 
The effective separation of oil-in-water dispersions involves high flux of water through the membrane 
and, at the same time, high rejection rate of the oil phase. The effects of transmembrane pressure 
and crossflow velocity on rejection of oil droplets and thin oil films by pores of different cross-section 
are investigated numerically by solving the Navier-Stokes equation. We found that in the absence 
of crossflow, the critical transmembrane pressure, which is required for the oil droplet entry into a 
circular pore of a given surface hydrophilicity, agrees well with analytical predictions based on the 
Young-Laplace equation. With increasing crossflow velocity, the shape of the oil droplet is strongly 
deformed near the pore entrance and the critical pressure of permeation increases. We determined 
numerically the phase diagram for the droplet rejection, permeation, and breakup depending of the 
transmembrane pressure and shear rate. Finally, an analytical expression for the critical pressure 
in terms of geometric parameters of the pore cross-section is validated via numerical simulations 
for a continuous oil film on elliptical and rectangular pores. 
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1 



I. INTRODUCTION 



With the recent advances in environmental and biological technologies, there has been in- 
creasing interest in characterization and modeling flows at the micron scales including flows 



in microchannels and nanochannels 



porous media 



Bfl, 



droplet-based microfluidics p, IsJ. The industrial app lications include oil extraction from 



HQ, 



multiphase flows through porous media |3|, 



and 



treatment of oily wastewater 



mm, 



and encapsulation of molecules, 



cells, and microorganisms (l3l-ll5|. In most of these processes, one or more phases are 



dispersed in a continuous phase in the form of emulsions, which are usually produced by 
shearing two immiscible phases against each other in the presence of surfactants [16]. In 
some cases, emulsions serve as means of transport of molecules, bio-reagents, and drugs, 
and ultimately provide the environment for enhanced reactions • Another techno- 

logical application of emulsions is to improve the transportability or displacement of highly 
viscous liquids. For example, heavy crude oil is emulsified to form a less viscous mixture 



to facilitate its transportation H8|, H9H. In addition, oil-in- water emulsions are used to en- 
hance recovery and increase sweep efficiency from crude oil reservoirs by blocking highly 
permeable paths and preventing channeling of the displacing fluid [9[. Common methods 
for separation of emulsions include evaporation of the continuous phase [20] . destruction 
(demulsification) [2]]], and membrane filtration flo| . 

Membrane microfiltration has proven to be an efficient way for separating oil-in-water 



emulsions 22 



4 l23[ . In comparison with the conventional methods of filtration (gravity sepa- 
rators, centrifuges, etc.), membrane microfiltration has several distinct advantages including 
reduced space requirements, higher permeate quality, and lower operating costs |10| . De- 
spite its advantages, microfiltration efficiency can be greatly reduced because of membrane 
fouling at highly concentrated emulsions or long filtration times j^j]. Fouling is generally 
caused by the accumulation of the rejected phase on the surface of the membrane or inside 
the pore. There are four main mechanisms (blocking laws) for membrane fouling, i.e., stan- 
dard blocking, complete blocking, intermediate blocking, and cake formation 25]. Complete 
blocking is common for very dilute mixtures and during the initial stages of filtration when 
some pores are sealed by droplets and particles, thus reducing the permeate flux |2a, \27 \. 
Accumulation of the rejected droplets on the membrane surface results in the formation of 
the so-called cake layer, which is sometimes referred to as the secondary membrane as it adds 
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3. 



This mechanism is dominant 



a hydraulic resistance to the microfiltration process [28- 
at the final stages of filtration when the water flux depends mainly on the thickness of the 
cake layer. 

The efficiency of the microfiltration process is determined by the properties of the mem- 
brane material and oil-in- water mixtures. For example, the permeate flux is highly depen- 
dent on the oil concentration, stability of the oil phase in water, and the size distribution of 
oil droplets fiol. I31I]. Moreover, the membrane properties such as membrane materia". 



10, 



, pore 



32(]. It 



size and morphology, and membrane geometry affect the permeate flow resistance 
was shown that slotted (rectangular) pores resulted in higher flux rates compared to circular 
pores for similar operating conditions because of the lower fowling rate of the slotted-pore 



membranes [33] . Another approach to reduce the fouling rate is to introduce crossflow above 
the membrane surface. This method, known as the "crossflow microfiltration" , reduces foul- 
ing by sweeping away the deposited drops and particles and, hence, decreases the thickness 
of the cake layer. Therefore, crossflow microfiltration systems tend to produce hig her per- 
meate fluxes for longer times compared to dead-end microfiltration systems |34j, [35|] . One of 
the aims of the present study is to investigate numerically the entry dynamics of oil droplets 
into a membrane pore in the presence of crossflow. 

The dynamics of droplet breakup in steady shear flow is determined by the relative com- 



petition of the viscous stress, pressure, and interfacial tension 



. In general, the breakup 



process is initiated by the droplet deformation, which is linearly proportional to the rate of 



shear 



371 ]. When the critical deformation is reached, the droplet assumes an unstable con- 



figuration and undergoes a transient elongation before it breaks up |36|] . It was also shown 
that the geometric confinement as well as the viscosity ratio of the dispersed and continuous 



phases influence droplet breakup 



38(| . In recent years, the problem of droplet deformation 



and breakup has 
boundary integral 



3een extensively studied numerically using Lattice Boltzmann 



43j, and Volume of Fluid (VOF) (44 



4611 methods. The VOF method 



used in the present study has proven to be a powerful and efficient interface tracking algo- 
rithm that is both conceptually simple and relatively accurate j47|] . Due to the conservative 
discretization of the governing equations in the VOF method, the mass of each fluid is ac- 
curately conserved 3- Furthermore, the ability of the VOF method to automatically 
capture local and global changes of the interface topology, e.g., coalescence and breakup of 
droplets, has made it attractive for various multiphase flow applications [471 ]. 
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During the last decade, a number of studies have investi gate d the process of droplet 



40. 
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49-51 . In these microfluidic 



formation using cross-flowing streams in T-shaped junctions 
systems, two immiscible liquids are driven through separate channels until their streams meet 
at a junction, where the dispersed liquid extends into the continuous stream, resulting in 
periodic formation of equal-sized droplets j^j]. Regardless of the specific channel geometry 
and wettability of the channel walls, breakup of the emerging droplet in a cross-flowing 
stream is determined by the viscous drag when the droplet remains unconfined by the 
microchannel [52|]. For unconfined T-junctions, it was demonstrated experimentally that 
the droplet size strongly depends of the crossflow rate of the continuous phase and only 
weakly on the flow rate of the dispersed phase |49|, |5C|. It was also shown that, for given 
value of continuous phase flow rate, the size of oil droplets decreases with increasing viscosity 
ratio of the oil and water jso| . 

In this paper, numerical simulations based on the Volume of Fluid method are performed 
to study the influence of transmembrane pressure and crossflow velocity on the entry dy- 
namics of thin oil films and droplets into pores of various cross-sections. We find that the 
formula derived in Ref. [53] for the critical pressure of permeation of an oil droplet into a 
circular pore agrees well with the results of numerical simulations. The numerical analysis is 
then extended to thin oil films covering pores with elliptical and rectangular cross-section in 
the absence of crossflow. In the presence of crossflow, we obtain numerically the phase dia- 
gram for the droplet rejection, permeation, and breakup as a function of the transmembrane 
pressure and shear rate, and study the details of the processes in three different regions of the 
phase diagram. These results are relevant to microfiltration of dilute oil-in-water emulsions 
at early stages before the formation of the cake layer. 

The rest of the paper is organized as follows. The details of numerical simulations are 
described in the next section. The analytical predictions based on the Young-Laplace equa- 
tion are reviewed in Sec. 1111 Al and verified numerically for an oil droplet on a circular pore 
in Sec. IIIIBi The critical pressure of permeation for pores with elliptical and rectangular 
cross-section is reported in Sec. IHICi The results for the oil droplet dynamics near circular 
pores in the presence of crossflow are presented in Sec. IIIIDl The conclusions are given in 
the last section. 
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II. DETAILS OF NUMERICAL SIMULATIONS 



Numerical simulations were carried out using the commercial software FLUENT [54]. 
In order to control the transmembrane pressure and the crossflow velocity, a user-defined 
function was written and compiled along with the main solver. The Volume of Fluid method 
was used to solve the multiphase flow problem [55]. For a two-phase fluid, this method is 
based on the fact that the two phases form an impenetrable interface, i.e., each cell is filled 
with either one of the phases (denoting a specific phase zone) or a combination of two phases 
(denoting the interface). This is achieved by introducing a variable a ) known as the "volume 
fraction" , which is defined as the ratio of the volume of fluid in the cell and the total cell 
volume; and it varies between and 1 



56, 



571 ]. An example of how the volume fraction 



varies near the interface is illustrated schematically in Fig.[TJ 

The interface is tracked by solving the transport equation for the volume fraction as 
follows: 

^ + V- (aV) = 0, (1) 

where V is the velocity vector. Equation (pp) states that the substantial derivative of the 
volume fraction is zero, and, therefore, the interface is convected by the velocity fields at 
the interface. After solving Eq. (pp), the material properties are computed by considering the 
fraction of each component in the cell; e.g., the density is given by 

p = ap 2 + (1 - a)pi, (2) 

where p is the volume-fraction-averaged density. 

One momentum equation is solved and the velocity field is shared between two phases as 
follows: 

-(pV) + V • (pVV) = -Vp + V • KVV + VV T )] + pg + F, (3) 

where g is the vector of gravitational acceleration, and F is the source term. In multiphase 
flow applications, the source term is the surface tension force per unit volume and it is 
non-zero only at the interface. Using the divergence theorem, the surface tension force is 
defined as the volume force in a cell as follows: 

F.^f^L. (4) 
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where a is the surface tension between two phases and k is the mean curvature of the interface 
in the cell. This force is related to the pressure jump across the interface (determined by the 
Young-Laplace equation) and it acts in the direction normal to the interface. The surface 
tension term tends to smooth out regions with large interface curvature [58| . If the interface 
is in contact with the wall, the normal vector (Va), which defines the orientation of the 
interface in the cell adjacent to the wall, is determined by the contact angle. The effect of 
static contact angle is taken into account by imposing the interface unit normal for a point 
(a cell in the Finite Volume method), n^, on the wall containing the interface as follows: 

= n w co8(9 st ) + n t sin(0 5t ), (5) 

where n w is the unit vector normal to the wall, is a vector on the wall and normal to the 



contact line, and 6 st is the static contact angle [59]. 

A SIMPLE algorithm was used for the pressure-velocity decoupling. The momentum 
equation was discretized using a second order upwind scheme. To reconstruct the interface 
and, consequently, solve the volume fraction transport equation, a PLIC (Piecewise Linear 
Interface Reconstruction) method was used [g^]. The pressure equation is discretized using 
a staggered mesh with central differencing. In FLUENT, the interfacial tension is modeled 



using the well-known model of Continuum Surface Force (CSF) of Brackbill et al [59f] . Using 
CSF, the surface tension volume force Eq. <^ is added as a source term to the momentum 
equation and the curvature is computed in terms of the vector normal to the interface n via: 

1 



Inl 



V)|n|-(V- n) . (6) 
n ; 



Interfacial effects in multiphase flows are described by the Young-Laplace equation, which 
relates the pressure jump across the interface to its mean curvature and the surface tension 
coefficient. For flows at the micron length scales, the viscous effects are dominant and the 
inertial effects are typically negligible. The capillary number is a measure of how viscous 
shear stresses are compared to the interfacial stresses and it is defined Ca = jiU /a ) where U 
is the characteristic velocity, \i is the fluid viscosity, and a is the surface tension coefficient. 

In the present study, the numerical simulations are performed to investigate a separation 
process of two commonly used liquids, i.e., kerosene and water. The density of kerosene 
is p = 889kg/m 3 and the viscosity ratio of kerosene and water at standard conditions is 
l^o/l^w = 2.4. It is assumed that the water is deionized; and, thus, the surface tension 
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coefficient a = 0.0191 N/m is used throughout the study [611 ]. Furthermore, we consider 



hydrophilic surfaces (for example poly vinyl- pyrrolidone |6l|]) with contact angles of kerosene 
in water greater than 90°. 

In our simulations, the mesh was generated in GAMBIT using the Cooper mesh scheme. 
This method works by sweeping the node patterns of specified source faces through the 
whole volume and the resultant mesh consists of an array of tetrahedral grids. For the results 
reported in the current study, we used about 30 cells along the pore diameter. To test the 
grid-resolution dependence, we considered 3 finer meshes that contained 50, 70, and 90 cells 
along the pore diameter. In the absence of crossflow, the simulations were performed for an 
oil droplet (r^ = 1.0 /xm) on a circular pore (r p = 0.2 /xm) at two transmembrane pressures 
(1.000 and 0.951 bar) slightly above and below the exact value of the permeation pressure 
0.976 bar predicted by the Young-Laplace analysis. In all cases, the droplet would either 
penetrate into the pore or reside at the pore entrance for at least 40 /xs. Furthermore, it 
was previously shown that the velocity of the contact line in the VOF method is inversely 
proportional to the logarithm of the mesh size [621], and, therefore, it is expected that 
the droplet velocity in the shear flow will depend on the grid resolution. However, in the 
present study, the oil droplet becomes temporarily pinned at the pore entrance by the 
transmembrane pressure, and thus the contact line velocity becomes much smaller than the 
flow velocity in the channel. Nevertheless, we have performed numerical simulations in the 
permeation, rejection, and breakup regions of the phase diagram with 4 times finer meshes 
and found that our results remain unchanged. 



III. RESULTS 



A. The Young-Laplace analysis for circular pores 

Effective separation of oil-in- water emulsions is controlled by several key parameters such 
as the membrane pore size, surface energy, size of oil droplets, surface tension, and pressure 
difference across the membrane. If the transmembrane pressure is relatively high, then oil 
droplets will most likely penetrate the membrane surface resulting in low rejection rates of 
the oil phase. On the other hand, low transmembrane pressures tend to limit flux of water 
through the membrane. Hence, the optimum operating conditions strongly depend on the 
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critical transmembrane pressure required for an oil droplet entry into a membrane pore. 

When the transmembrane pressure across a hydrophilic membrane exceeds a certain criti- 
cal value, the oil phase will penetrate the membrane. Thus, for high separation efficiency, the 
transmembrane pressure should be maintained at a value below P cr iu which for a continuous 
oil film on the membrane surface with circular pores is given by 

_ 2acosfl 

r crit -> v ' / 

where a is the surface tension coefficient between oil and water, 9 is the contact angle of the 
interface of an oil droplet on a membrane surface immersed in water, and r p is the membrane 



pore radius 



531 ]. The critical pressure in Eq. ([7j) is determined by the Young-Laplace pressure 



due to the curvature of the oil-water interface inside the pore. 

If instead of a thin oil film, a droplet of oil is placed at the entrance of the membrane pore, 
the formula for the critical pressure, Eq. ([7j), has to be corrected by a factor that depends on 
the ratio Vd/r p to include the effect of the oil- water interface curvature above the membrane 
surface. It was previously shown [53|, [6jJ| that the pressure required to force an entry of an 
oil droplet of radius into a circular pore is given by 



_2<jcos# I 2 + 3cos#-cos 3 # 

cnt ~ rp y 1 ~ 4 (r d /r p f cos 3 # - (2 - 3 sin# + sin 3 #) ' [&) 



In the limit — >> oc, the curvature of the droplet above the pore vanishes, and thus this 
formula corresponds to a continuous oil film on the membrane surface, and Eq. (JSj) converges 
to Eq. ([71). Contrary to the case of the thin film, the critical pressure for an oil droplet with 
9 = 90° is negative, and the droplet will penetrate into the pore in the absence of the 
applied pressure gradient. We also comment that no such expression exists for the case 
when crossflow is present and the hydrodynamic drag force is exerted on the droplet parallel 
to the membrane surface. 

In what follows, we investigate the dynamics of an oil droplet and thin film entry into 
pores with various cross-sections. The critical pressures [Eq. ([7]) and Eq. (jSJ)] are compared 
with the results of numerical simulations in order to validate the numerical scheme. The 
numerical analysis is then extended to pores with rectangular and elliptical cross-sections. 
Finally, the effect of shear flow on the droplet entry is considered and a phase diagram of 
the transmembrane pressure versus shear rate is determined numerically. 
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B. An oil droplet on a circular pore in the absence of crossflow 

We first present the numerical results for the critical pressure required to force an entry 
of an oil droplet into a cylindrical pore. In the numerical scheme, the oil droplet is initially 
immersed in water above the membrane surface in the absence of flow. The transmembrane 
pressure is then gradually increased to a predetermined value. As the simulation continues, 
the droplet approaches the membrane surface and resides at the pore entrance. Depending 
on the applied pressure difference across the membrane, the droplet will either remain at the 
pore entrance or penetrate into the pore. We comment that when the applied pressure is 
close to the critical pressure, the dynamics of an oil droplet entry into the pore is significantly 
slowed down because the net driving force on the droplet is reduced. 

The critical pressure as a function of the droplet radius is plotted in Fig. [2] using Eq. (JSj) 
for three values of the pore radius. The error bars in Fig. [2] indicate the upper and lower 
values of the transmembrane pressure when the oil droplet either enters the pore or remains 
at the pore entrance during the time interval of about 40 /xs. We find an excellent agreement 
between the results of numerical simulations and analytical predictions of Eq. (jSJ) , which 
provides a validation of the numerical method. Furthermore, as shown in Fig. El the critical 
pressure increases with increasing droplet radius or decreasing pore radius. As the droplet 
radius increases, the critical pressure approaches an asymptotic value predicted by Eq. ([Tj) 
for a thin oil film (not shown) . 

C. Critical pressure for pores with arbitrary cross-section 

We next investigate the influence of pore cross-sectional shape on the critical transmem- 
brane pressure using simple physical arguments and numerical simulations. According to 
the Young-Laplace equation, the pressure jump across an interface between two immiscible 
fluids is related to its mean curvature and the surface tension as follows: 

AP = 2(jk, (9) 

where k is the mean curvature of the interface computed by averaging two principle curva- 
tures. In the absence of gravity, the mean curvature of an arbitrary surface z{x, y) is given 
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by 

2* = V-(-=EiL=Y (10) 

It follows from Eq. (|9j) that an interface, which is subject to a prescribed pressure jump and 
constant surface tension coefficient, has a constant mean curvature. Therefore, if the gra vity 



is negligible, the fluid-fluid interface forms a section of the so-called Delaunay surface |63l. 164] . 

As shown in Fig. El if the oil-water interface is bounded by the walls of a pore of arbitrary 
cross-section, the constant contact angle at the pore surface imposes a boundary condition 
for Eq. f [T0l ) in the form 

cos<? = n-(-= (11) 



Vi + |v^| 2 

where n is the outward unit vector normal to the pore surface [6. 

Integrating Eq. ( fTUl) over an arbitrary cross-section with a smooth boundary and using 
the divergence theorem, we obtain 

V2«= / n- ( — V l = )dL, (12) 



!c P \y/i + \Vz\ 2 < 

where C p and A p are the perimeter and cross-sectional area respectively. Taking the integral 
on the right hand side over the perimeter gives the following relation for the mean curvature 
of the oil- water interface |66| 

2.= ^-. (13) 

According to Eq. f fT3l ) , the mean curvature of an interface bounded by a pore of arbitrary 
cross-section can be related to the geometric properties of the boundary. Remember that 
Eq. (|9j) relates the mean curvature to the surface tension and the pressure jump. In the case 
of an oil film on a pore with an arbitrary cross-section, the critical applied pressure is equal 
to the pressure jump at the interface. Therefore, combining Eqs. (|9j) and ( TT31) . we obtain the 
critical permeation pressure for the oil film to enter into a pore of arbitrary cross-section 

Pm = tSpl. (14) 

This equation can also be derived from the force balance between the applied pressure and 
the Laplace pressure due to the curvature of the oil-water interface inside the pore. 

It should be noted that the boundary value problem Eq. f fTUl ) does not always have a solu- 
tion for a stable interface with a constant mean curvature and a constant contact angle 
In other words, there is a limitation on the values of the contact angle that correspond to 
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the attached interface with a constant mean curvature. For example, if the contact angle 
(computed from the oil phase) is larger than the critical value, then Eq. (fTOl) subject to 
the boundary condition Eq. ( TTTT) does not have a stable solution. As a result, the interface 
cannot remain attached to the bounding surface with a prescribed contact angle and, at 
the same time, maintain a constant mean curvature required by the Laplace equation. The 
critical contact angle is determined by the largest curvature of the cross-sectional shape for 



smooth boundaries and by the smallest wedge angle for boundaries with sharp corners [681 ]. 

In what follows, we consider two special cases of rectangular and elliptical pores and 
compare predictions of Eq. ( TH1) with the results of numerical simulations. The problem is 
illustrated schematically in the inset of Fig.[H The oil film covers the pore entrance of a 
hydrophilic membrane subject to a pressure gradient (the transmembrane pressure). Since 
the Bond number is small, Bo = (p w — p Q ) r^g/ a = 6 x 10 -8 , the effect of the gravitational 
force can be neglected. 

1. Thin oil film on the membrane surface with a rectangular pore 

In this subsection, we investigate the dynamics of an oil film entry into a rectangular 
pore, which is sometimes referred to as a "slotted pore" [69]. The rectangular shape of a 
slotted pore is characterized by the width (the shorter side) and the length (the longer side) . 
Using A p = vol and C p = 2(w + /), the corresponding critical pressure is obtained from 
Eq. (pi]) as follows: 

P cmt = 2acos0(- + \), (15) 



<w I J' 

where w and I are the width and length of the pore cross-section [66|, [70, IZll] . In the limit 
when / 3> Eq. f |T5l ) reduces to P crit = 2 a cos 9 / w, which is the critical pressure on an oil 
film entering into an infinitely long rectangular pore. In this case, one of the curvatures of 
the interface is zero and the other curvature is proportional to the width of the pore, and 
thus the shape of the interface is a part of a cylinder with the radius w/2cos9. 

The results of numerical simulations and predictions of Eq. ( |T5l) are shown in Fig. [4] for 
several aspect ratios. As expected, square pores have the highest critical pressure due to 
the largest perimeter-to-area ratio. The critical pressure decreases with increasing aspect 
ratio. These results demonstrate that there is an excellent agreement between the numerical 
results and analytical predictions based on the Young-Laplace equation. It is important to 
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note that when the applied pressure is close to the critical pressure, the net force on the 
interface is small, and, therefore, very long simulation time is required to capture the motion 
of the interface. The symbols in Fig. 0] indicate the rejection and permeation pressures that 
were resolved numerically without excessive computational effort. Interestingly, each curve 
in Fig. H] is well described by the function P crit = 2 a cos 9 / /, which is shifted upward by a 
constant 2 a cos 6 / w (indicated by the horizontal lines in Fig.Hj). 

It was previously shown that for a rectangular cross-section of arbitrary aspect ratio, 
Eq. ( ITUl) has a solution with a constant curvature for a non- wetting fluid (9 > 90°) when 



the the contact angle 9 < 135° 



68, 



71 



l73l| . In the case of a square pore, the interface is part 
of a sphere with the radius w/2cos For other aspect ratios, the interface surface 

has a constant mean curvature k = cos# + but it is no longer spherical fssl ] . 

Figured shows snapshots of the oil- water interface inside the square pore obtained from 
our numerical simulations. The transmembrane pressure is set to a value computed from 
Eq. (PSD for the contact angles 9 = 90°, 120°, 135°, and 150°. As observed in Fig. El the 
concave shape of the interface strongly depends on the contact angle. When 9 = 90°, the 
interface enters the pore with zero curvature and, according to Eq. f fT5l) . with zero pressure 
gradient. For contact angles between 90° and 135°, the interface bends in the center and 
penetrates into the pore before it starts to move near the corners. For 9 > 135°, the distance 
between the interface location in the center and at the corners will theoretically be infinity 
because the interface becomes pinned at the corners while the inner part penetrates into the 
pore j75l ]. 



2. Thin oil film on the membrane surface with an elliptical pore 

The critical permeation pressure for an oil film covering a pore with an elliptical cross- 
section can be estimated from Eq. ( IT41) and the geometric properties of an ellipse. However, 
there is no exact expression for the perimeter of an ellipse. In our study, we use one of the 
most accurate and compact approximations that predicts the perimeter of an ellipse with 
an error of -0.04% Q 

Q 7 

C p ^n(a + b) 1 + , (16) 

P V 7 L 10 + V4-3/J 
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where a and b are the major and minor radii of the ellipse and h — (a — b) 2 /{a + b) 2 . Using 
Eq. (|T6l) and the expression for the ellipse area A p = n a 6, the critical pressure is given by 



z> ~ + 

crit ~ 7 

ab 



3h 

1 + 



acos#. (17) 



10 + V4-3/J 

Clearly, in the case of a circular pore, a = b = r p , Eq. (JTTJ) is reduced to Eq. ([7j). 

The results of numerical simulations and predictions of Eq. ( 1T71) are summarized in Fig. [6] 
for different aspect ratios. Similar to the case of the rectangular pore, the symbols indicate 
pressures of rejection and permeation that were resolved during the simulation time interval 
of about 40 lis. It can be observed that the critical pressure decreases with increasing ellipse 
aspect ratio, and the numerical results agree very well with predictions of Eq. ( fTTl) . We also 
comment that when a>6, each curve in Fig. [6] asymptotes to P crit w 4 a cos 9/n 6, which is 
higher than the value P crit = 2 a cos / w estimated for an infinitely long rectangular pore 
(see section MIC 11) . This difference arises because an infinitely long ellipse and an infinitely 
long rectangle of equal width have the same perimeter but different areas. 

If the aspect ratio of an elliptical pore is less than 1.635, then the boundary value problem 
given by Eq. ( JTU1) has a solution for any contact angle jizj]. However, when a/b > 1.635, there 
is a critical contact angle above which Eq. ( ITOl) has no solution [77| . For a/b = 1.635, the 
critical contact angle is 180°, and it decreases with increasing aspect ratio. The largest 
aspect ratio considered in the present study is a/b = 4.0 for which the critical contact angle 
is 153.05° [77j. Therefore, the contact angle of 120° used in our simulations generated an 
interface with a constant mean curvature described by Eq. ( |T3l ) . We comment that the oil- 
water interface inside the elliptical pore is not spherical because the boundary condition 
Eq. f lTTl) is not satisfied at the intersection of a sphere and a cylinder with an elliptical cross- 
section. Similar to the rectangular cross-section, we found that at the critical pressure given 
by Eq. f lTTl) . the contact line is pinned at the antipodal points of the highest curvature of 
the ellipse, (x, y) = (±a, 0), while the rest of the interface penetrates into the pore. 



D. Sheared droplet on the membrane surface with a circular pore 

We next examine the combined effect of the transmembrane pressure and crossflow ve- 
locity on the entry dynamics of an oil droplet into a circular pore. The computational setup 
is illustrated schematically in Fig. [71 The shear flow is induced by translating the upper 
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wall with a constant velocity. In our simulations, the effective shear rate 7 is defined as 
the ratio of the upper wall velocity to the channel height. The relevant dimensionless num- 
bers, the capillary and Reynolds numbers, are estimated to be Ca = iiw^r^ja < 0.03 and 
Re = Pw^d/l^w ^ 0.5. The width of the channel is chosen to be about 8 times larger than 
the droplet radius in order to minimize finite size effects in the lateral direction. Initially, 
the droplet is released upstream to insure that the flow reaches a steady state before the 
droplet approaches the pore. At the same time, the transmembrane pressure is set to a 
prescribed value, and the simulation continues until the droplet either reaches the outlet, 
penetrates into the pore, or breaks up. 

The main results of this study are summarized in Fig. [HI which shows the phase diagram 
for the droplet rejection, permeation, and breakup depending on the transmembrane pressure 
and shear rate. The corresponding snapshots of the droplet for five different cases (denoted 
by the capital letters A, B, C, D, and E) are presented in Fig.[9l Below, we discuss the 
details of the processes in the three different regions of the phase diagram and provide an 
estimate of the leakage volume during the droplet breakup. 

In the permeation region shown in Fig. El the transmembrane pressure is larger than the 
steamwise drag force, and the oil droplet penetrates into the pore. A series of snapshots at 
point B in Fig.[9]demonstrate the details of the permeation process. As observed in Fig.[9](B), 
at first, the droplet partially penetrates into the pore and becomes strongly deformed in the 
shear flow. However, the droplet does not breakup because its size above the pore decreases 
as the droplet penetrates into the pore, and the viscous shear stress acts on a progressively 
smaller surface area. 

With increasing shear rate, the effect of viscous forces becomes more important, resulting 
in strong deformation of the droplet shape near the pore entrance. We find that at sufficiently 
large transmembrane pressures and Ca > 0.015, the oil droplet breaks up. In this case, the 
larger droplet is washed off downstream and the smaller droplet enters the pore, and, as a 
result, the membrane leaks. 

Depending on the values of the transmembrane pressure and shear rate, two different 
breakup regimes were observed. The first regime is bounded by the minimum breakup 
pressure, which is found from Fig. [8] to be 1.00 ± 0.01 bar. Above this pressure, a small 
fragment is detached and penetrates into the pore while the main droplet is carried away by 
the shear flow. During this process, the droplet has a limited time to deform, and, therefore, 
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the breakup occurs only locally without significant deformation, as shown in Fig.[9](D). 

The breakup process at higher transmembrane pressures occurs in a qualitatively different 
way [see Fig. 0(E)] . As the droplet approaches the pore, it momentarily slows down and 
remains at the pore entrance during the "residence time". In this case, the effects of the 
drag and the transmembrane pressure are relatively large and comparable with each other. 
As a result, the shape of the droplet is significantly deformed by the shear flow and a thin 
bridge is formed between two parts of the stretched droplet. This thinning is known as 



m 



791 ] . For a short 



necking and it usually indicates the initial stage of the breakup process 
time interval, the thin neck holds the two parts of the droplet together. At the final stage 
of breakup, the neck gets thinner and thinner near the the edge of the pore, and at some 
point, it becomes unstable and the droplet breaks. 

Visual inspection of the snapshots of the droplet near the pore entrance revealed that, at 
a given shear rate in the breakup regime in Fig. El the residence time is roughly independent 
of the transmembrane pressure. Therefore, it is expected that during the breakup process, 
the volume of leaked droplets is proportional to the applied pressure. Figured!)] shows the 
leakage volume as a function of the transmembrane pressure when 7 = 5 x 10 5 s _1 . Indeed, 
the leakage volume is almost linearly proportional to the applied pressure, indicating that 
the flow inside the pore is described by the Hagen-Poiseuille equation. 

The lower part of the phase diagram in Fig. [8] indicates operating conditions when the oil 
droplet is rejected by the membrane and washed off by the shear flow. An example of the 
rejection process at low shear rates is presented in Fig.[9](A). Although the droplet partially 
penetrates into the pore, the flow generates a force on the droplet surface, which pulls the 
droplet away from the pore, resulting in the droplet rejection. As shown in Fig.[9](C), the 
residence time at higher shear rates is reduced, and the droplet is carried away by the flow 
without penetrating into the pore. 

The threshold of permeation is determined by the competition between the drag force and 
the transmembrane pressure. Naturally, with increasing shear rate, the drag force increases; 
and, therefore, it is not surprising that the boundary curve separating the permeation and 
rejection regions in Fig. [8] increases with shear rate. However, it is difficult to estimate 
the exact dependence of the drag force on the droplet because its shape becomes strongly 
deformed in the shear flow. We also comment that, in the range of shear rates reported in 
Fig. [HI the lift force is about an order of magnitude smaller than the drag force [49 ]. 
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The critical shear rate that marks the boundary between the permeation and breakup 
regions in Fig. [8] can be estimated using simple force balance arguments. In the absence of 
gravity, the sheared droplet is subject to forces of surface tension, Laplace pressure, drag, 
and lift. Neglecting the lift force and following the analysis in Ref. j49|, the torque balance 
around the edge of the pore for the droplet configuration depicted in Fig.[9](E£ 6 ) can be 
written as follows: 

F D d d + (P 2 -P 1 )A n d n -F a d n = 0, (18) 

where d n and A n = nd^/A are the diameter and cross-sectional area of the thinnest stable 
neck, F a = tt a d n is the surface tension force around the perimeter of the neck, and P 2 and 
Pi are the pressures inside the droplet and in the channel, respectively. The Stokes drag 
force on the spherical droplet, Fp 1.1 n ji w 70! J, is estimated for the viscosity ratio 2.4 and 
the average flow velocity 7^/2 [50|]. In our simulations, the typical diameter of the thinnest 
stable neck is d n « 0.9rf p [see Fig.[9](E£ 6 )]. Using d d = 1.7/xm, a = 0.0191 N/m, and \i w = 
10 _3 kg/ms in Eq. (fl8l) . the critical shear rate is roughly estimated to be 7 ~ 3.4 x 10 5 s _1 , 
which is in good agreement with the value 7 w 3.2 x 10 5 s _1 obtained numerically in Fig. [51 
The permeation, rejection, and breakup regions identified in the phase diagram in Fig. [8] 
can be useful for the optimal design and operation of crossflow microfilt ration systems. It is 
apparent that the permeation region should be avoided for filtration purposes. The optimal 
performance of the microfiltration system with maximum rejection is achieved in the upper 
part of the rejection region where the large transmembrane pressure results in high flux of 
water while the oil phase is completely rejected. However, the separation efficiency can be 
increased at higher transmembrane pressures in the breakup region, where the higher flux 
of water is accompanied by some oil leakage. 



IV. CONCLUSIONS 



In this paper, the numerical simulations were carried out to investigate the influence of 
the transmembrane pressure and crossflow velocity on the entry dynamics of thin oil films 
and oil droplets into pores of different cross-section. We considered hydrophilic membrane 
surfaces with contact angles of oil in water greater than 90°. The numerical method was 
validated against the analytical solution for the critical pressure of permeation of an oil 
droplet into a circular pore in the absence of crossflow. Furthermore, we found that the 
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results of numerical simulations of thin oil films on elliptical or rectangular pores agree 
well with the theoretical prediction for the critical pressure expressed in terms of geometric 
parameters of the pore cross-section. Also, examples of curved oil-water interfaces inside 
elliptical and rectangular pores were discussed for different aspect ratios and contact angles. 

In the presence of crossflow above the membrane surface, we have determined numerically 
the phase diagram for the droplet rejection, permeation, and breakup as a function of the 
transmembrane pressure and shear rate. A detailed analysis of the droplet dynamics near the 
pore entrance was performed in the three different regions of the phase diagram. We found 
that in the permeation region, the transmembrane pressure is larger than the steamwise drag, 
and the oil droplet penetrates into the circular pore. With increasing crossflow velocity, the 
shape of the droplet becomes strongly deformed near the pore entrance; and, at sufficiently 
high transmembrane pressures and shear rates, the droplet breaks up into two fragments, 
one of which penetrates into the pore. It was also shown that during the breakup process, 
the residence time of the oil droplet at the pore entrance is roughly independent of the 
transmembrane pressure, and the volume of the leaked fragment is nearly proportional to 
the applied pressure. Finally, the numerical value of the critical shear rate that separates 
the permeation and breakup regions is in good agreement with an estimate based on the 
force balance arguments. 
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FIG. 1: (Color online) Values of the volume fraction for each of the phases and the interface (a = 1 
in the first phase, a — in the second phase, and < a < 1 at the interface). 




FIG. 2: (Color online) The critical pressure computed using Eq. (|8j) for the values of the pore 
radius 0.15 /im (continuous line), 0.2 /xm (dashed line), and 0.3 /im (dotted line). Error bars are 
extracted from the numerical simulations (see text for details) with the parameters p /Pw — 2.4, 
Po/pw — 0.781, a = 0.0191 N/m, and 9 — 135°. The symbols indicate (x) rejection and (o) 
permeation of the oil droplet. 
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FIG. 3: Schematic representation of the oil-water interface inside the pore of arbitrary cross-section 
with perimeter C and area A. The interface forms a constant angle 9 with the inner surface of the 
pore. 
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FIG. 4: The critical permeation pressure for the oil film into the rectangular pore with different 
aspect ratios. The curves are Eq. ([T5]) and the symbols are the numerical results for Po/Pw — 2.4, 
Po/pw = 0.781, a = 0.0191 N/m, and 9 = 120°. The symbols indicate (x) rejection and (o) 
permeation of the oil film. 



25 



(a) = 90° (b) = 120 




(c) = 135° (d) = 150° 

FIG. 5: (Color online) Snapshots of the oil- water interface inside the square pore for contact 
angles 9 = 90°, 120°, 135°, and 150°. The critical pressure Eq. (TT5]) is computed for the surface 
tension a = 0.0191 N/m and the pore width 1.5 fim. Cases (a), (b), and (c) correspond to the 
stationary interface, while in the case (d) the interface is in transient state (see text for details). 
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FIG. 6: The critical pressure of permeation for the oil film into the elliptical pore as a function 
of the major and minor radii. The curves are computed using Eq. (jT7|) . The symbols represent 
numerical results for the parameters fi / f fi w = 2.4, p /p w — 0.781, a = 0.0191 N/m, and 9 — 120°. 
The symbols indicate (x) rejection and (o) permeation of the oil film. 
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Pressure Outlet 



FIG. 7: (Color online) Schematic representation of the oil droplet in the channel with the circular 
pore. The shear flow is induced by the upper wall moving with a constant velocity parallel to 
the stationary lower wall. The droplet radius is = 0.9 fim and the pore radius is r p = 0.2 fim. 
Periodic boundary conditions are applied at the inlet and outlet, while a constant pressure is 
maintained at the side walls. 
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FIG. 8: (Color online) The phase diagram of the transmembrane pressure versus shear rate for the 
oil droplet with = 0.9 fim and the circular pore with r p = 0.2 /im. The contact angle is 9 = 135°. 
Each symbol represents a separate simulation that corresponds to either (O) permeation, (□) 
breakup or (V) rejection. Letters A, B, C, D, and E indicate operating conditions for the series of 
snapshots shown in Fig.[9l 
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FIG. 9: (Color online) Sequences of snapshots of the sheared droplet on the pressurized pore 
for different operating conditions as indicated in Fig.[8l The droplet radius is — 0.9 /im, the 
pore radius is r p = 0.2 /im, and the contact angle is 9 = 135°. The letters denote (A) rejection 
at low shear rates (At « 15 /is), (B) permeation (At « 15 /is), (C) rejection at high shear rates 
(At « 2/is), (D) local breakup (At « 2.5 /is), and (E) breakup with necking (At « 4/is). 



30 




i Lj I I I I I I I I I I I I I I I I I I I 

0.05 0.1 0.15 0.2 

Volume of Leakage (|um 3 ) 



FIG. 10: (Color online) The leakage volume as a function of the applied pressure for = 0.9 /im, 
r p = 0.2 /im, and 7 = 5 x 10 5 s _1 . The square symbols indicate the numerical results and the 
dashed line is the best fit to the data. The inset shows a snapshot of the process shortly after 
breakup of the droplet. The contact angle of the oil droplet in water is 9 = 135°. 
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